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Abstract 

We propose a ferromagnetic/ferroelectric hybrid double quantum well structure, and present an 
investigation of the Curie temperature (T c ) modulation in this quantum structure. The combined 
effects of applied electric fields and spontaneous electric polarization are considered for a system 
that consists of a Mn 5-doped well, a barrier, and a p-type ferroelectric well. We calculate the 
change in the envelope functions of carriers at the lowest energy subband, resulting from applied 
electric fields and switching the dipole polarization. By reversing the depolarizing field, we can 
achieve two different ferromagnetic transition temperatures of the ferromagnetic quantum well in 
a fixed applied electric field. The Curie temperature strongly depends on the position of the Mn 
5-doped layer and the polarization strength of the ferroelectric well. 
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The diluted magnetic semiconductor (DMS) has been generally known as one of the 
promising candidates for spintronic device materials in virtue of the coexistence of ferro- 
magnetic and semiconducting properties in it. Ferroelectric material has also attracted 
significant interest because of its promising potential in various technological applications, 
such as binary data storage media in nonvolatile random access memories due to its spon- 
taneous electric polarization. In both research fields, many experimental and theoretical 
studies have been performed. 

Because the spintronic devices should ultimately be operated at room temperature, much 
effort has been focused on increasing the ferromagnetic transition temperature (T c ) of DMS 
above room temperature. Among many materials, ZnMnO is considered to have T c above 
300 K with 5% Mn per unit cell and 3 x 10 20 holes per cm 3 according to a theoretical 
prediction. [j| 

Recently, observation of the ferroelectric properties was reported in Li-doped ZnO bulk 
samples. 0, [| The reason for the ferroelectric property is attributed to the following: when 
the size of the dopant Li atom (0.6 A) is smaller than the host Zn atom (0.74 A)[2|, then 
the Li atoms can occupy off-center positions, thus locally inducing electric dipoles, thereby 
leading to ferroelectric behavior like PbGeTe ^| ferroelectric semiconductor (FES). 

The mechanism of DMS ferromagnetism is classified upon materials and growing tech- 
niques. The first class of approach is ferromagnetism due to the Ruderman-Kittel-Kasuya- 
Yoshida(RKKY)/Zener indirect exchange interaction by delocalized holes (hole mediated) 



based on the mean-field approximation, p, |f| The second class of approach is also carrier- 
induced ferromagnetism as a results of KKR-CPA-LDA (Korringa-Kohn-Rostoker coherent- 
potential approximation and local density approximation) calculations of the electronic 
structure of doped DMS alloys. [?], 0] The third class of approach suggests the hole hopping 
mediated ferromagnetism between polarons having strongly localized charge carriers. 3, |H| 
And the fourth one is ferromagnetism due to the ferromagnetic clusters or secondary 
phases. Therefore, it is necessary to decide on a case-by-case basis which mecha- 

nism is applicable. In our work, we apply the first class of ap pro ach based on the mean-field 



theory for carrier- induced ferromagnetism in a DMS. 
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Using ideas based on the dependence of T c of DMS on the spatial distribution of mag- 
netic ions, and envelope functions of carriers at the lowest energy subband in a confining 
potential|lfiL llfiL Il7j| . we model a hybrid double quantum well (HDQW) system shown in Fig. 
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1(a). The structure of ZnO/Zni_. E Mg.j;0/Zni_< y Li, / has the upper p-type ZnO well with an 
additional Mn 5-doped layer at the middle of the well (or at the upper edge of the well). 
The p-type might be achieved by the doping of group V 18, 2 or group ij^J elements. 
The lower well Zn^LiyO is the p-type ferroelectric well with spontaneous polarization P. 
The inverse potential profile of a hole is shown in Fig. 1(b), for dipole up and dipole down 
cases respectively. Because screening lengths at the interfaces between ZnLiO and ZnMgO 
are assumed to be different, the potential profile is asymmetric (different values of Vi and 
V2) in region W2. The dimension of the structure are chosen: Wi = 10 nm, W2 = 10 
nm, B=5 nm, and the capping layer is 20 nm. The confinement potential is Vo=-263 meV 
with 20% of Mg per unit cell in ZnMgO barriers. |2l| Because our work is applicable in a 
regime of low carrier density, occupying only the lowest energy subband of a heavy hole, 
carrier concentrations of both wells are in the order of 10 11 /cm 2 [l5| and the additional band 



bending due to carriers is small. 22] 



We previously demonstrated electric field control o 



ferromagnetism in Mn 5-doped asym- 
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but present structure is different 



metric conventional double quantum wells (CDQW) 
from CDQW because of hybridizing ferroelectricity. In this work, we can obtain additional 
effects from the ferroelectric well due to the reversal of spontaneous depolarizing fields. We 
can control the number of holes around the Mn 5-doped layer by reversing depolarizing fields 
as well as applying electric fields. The main purpose of this work is to show the possibility 
of using ferromagnetic/ferroelectric hybrid structures to modulate T c by reversing polariza- 
tion. We can obtain two different ferromagnetic transition temperatures in the same applied 
electric field in HDQW. 

The Hamiltonian for this system is given by 

H = -d z ^—d z + V c (z) - F g z + V d (z). (1) 
m{z) 

Throughout this calculation, we adopt atomic units R = m e 4 /2h 2 for the energy unit, and 
Ob = h 2 /moe 2 for the length unit, where tuq is the free electron mass. Here m! = m* /rriQ 
and m* is the hole effective mass and F g is the carrier charge times for the applied electric 
fields. Vd is the potential due to the spontaneous depolarizing field. Because we would like to 
provide a qualitative estimation without complication, we ignore the Hartree and exchange- 
correlation interactions among carriers, and self-consistent solving of electrostatic potential 
with the Poisson equation even though the Hartree and exchange-correlation interactions 



are related to the effective mass 



22j and distribution of carrier spins. (^J The thickness of 
the barrier is enough so that we can ignore coupling by tunneling between the wells. We 
attempt to solve the eigenvalue problem, Hip(z) = Eif)(z), and it becomes 

—d 2 Mz) - [V c {z) - F g z + V d (z) - E)iP(z) = 0. (2) 

The carrier confinement potential V c (z) is Vq inside the wells and outside the wells. The 
potential due to dipoles is Vd(z) = F d x [z — B/2) + (Vq ± Vi) where Fa is the spontaneous 
depolarizing field and Vi(> 0) is the electrostatic potential at the interface due to screening 
charges by the Thomas- Fermi model of screening. ± signs correspond to dipole left and 
dipole right cases. This potential profile is shown in Fig. 1(b) with F g = 0. We know that 
the general solution for Eq. (J2J) is a linear combination of Airy functions, Ai(x) and Bi(x). 
Thus, we write 

if>(z) = C x Ai{i) + C 2 Bi{i), (3) 

where 



i = ±{m'\F g + F i \ r l*[z- 



E - V =f V x + F t 
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9 '- u L ~ (F g + F d) 

with ± = Sgn(F g + F d ) for the region W 2 and 

for elsewhere with F g > assumed. By using boundary conditions and the continuity of 
wave functions and their derivatives at the boundaries of each region, we can calculate C\ 
and C*2 and energy eigenvalues in the system numerically. In principle, we should include 
the effects of the exchange interaction between Mn ions and carriers on the carrier wave 
function^ 0, but we ig no r e these effects beeause our SyS te m bas only a ve ry thin M„ 
5-doped layer. Therefore, our wave function is limited to a system with very thin Mn layers 



submonolayer) . 26| We write T c in the form 
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Tc = ^T^^ (4) 
LZkb nn, J 

Here, c(z) is a magnetic ion distribution function, J vd is the exchange integral of carrier-spin 
exchange interaction, and S is a Mn ion spin. We calculate the change in the fourth power 
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of growth direction envelope functions, ^(z)! 4 , of carriers at the lowest energy subband in 
the HDQW as a function of the applied electric fields. 

In numerical calculations, we choose physical parameters for p-type 
ZnO/Zn 1 _ a ,Mg :r O/Zn 1 _j / Lij / 0, with x = 0.2 and y = 0.05 The confinement potential 



Vq = —263 meV, and m' h = 0.78 and 1 for the well and the barrier, respectively. [21, Q In 
our calculation, the exact values of J p d and S are not required because we calculate the 
ratio of T c in Eq. (j3J). Usually J p d is one of the important factors in determining the size of 
T 

- 1 c- 

Figure 2(a) shows the dependence of the ratio of ferromagnetic transition temperatures, 
T c /T c0 , on bias voltage applied across the dipole left HDQW, for Mn center-doped (open 
triangle), and Mn edge-doped (closed circle) respectively. Here T c0 is T c at F g = meV/nm 
for the Mn center-doped CDQW without Li doping, [l^ For Zno.95Lio.05O, there are approx- 
imatelyl.05 dipoles/nm 3 , and they induce a maximum polarization of approximately 8.77 
fiC/cm 2 . While there is no depolarizing field (a full screening case) when barriers are metal, 
the maximum depolarizing field can be up to 1.2 eV/nm[2^| when barriers are insulators 
(no screening charge). Our case is in between those cases, the depolarizing field Fd = —0.01 
meV/nm and the screening electrostatic potential Vi=0.02 meV are used as input parame- 
ters for dipole down FES. Applied electric fields shift the envelope functions according to the 
change of potentials, as shown in Fig. 2(b). Then, the effective hole concentration increases 
(decreases) in the Mn edge-doped layer (center-doped), and T c increases (decreases) as a 
result. When F g is less than 2 meV, there are few carriers confined at the lowest subband 
of the Mn doped well because of the effect of negative (dipole left) Fd. 

Figure 3 shows the dependence of the ratio of ferromagnetic transition temperatures 
T c /T c0 on bias voltage applied across the Mn center-doped HDQW for dipole down (closed 
circle), and dipole up (open triangle) respectively. By reversing the direction of spontaneous 
polarization, the change in Curie temperature occurs below F g =2 meV. This effect is caused 
by asymmetry of electrostatic potential due to screening charges. The larger asymmetry is, 
the more effective the reversal is. Therefore, it is important to fabricate a sample having 
asymmetric potential to obtain this result. Both dipole down and dipole up cases have the 
same value of T c /T c0 above F g =2 meV, because the envelope functions of carriers depend 
on the potential profile of the Mn doped well (left side well), which are not affected by Fd 
as shown in the inset. When the coercive field of the ferroelectric well is much smaller than 
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the F g =2 meV, the Curie temperature may change at the coercive field. But we do not take 
the reversal of polarization due to F g into account, because we are interested in the regime 
of F g lower than the coercive field. The inset shows the change in the potential profile due 
to the reversal of dipole polarization in the HDQW as a function of z at F g =4 meV/nm. 
The potential profiles in the left-hand side well are the same regardless of the direction of 
dipoles. The energies of the lowest subbands are the same. These energy eigenvalues are 
shown in Fig. 4. 

Figure 4 shows energy eigenvalues as a function of applied electric fields for HDQW. Solid 
lines indicate the upper energy limit of the Mn delta-doped well. The energy degeneracy 
occurs at F g = for CDQW. When we compare the dipole right case and dipole left case 
of Fig. 4, the energy degeneracy at F g = shifts to the left (dipole right) or to the right 
(dipole left), because the energy levels corresponding to the Mn doped well are the same, 
and only the energy levels of the FES well are shifted up (dipole right) or down (dipole left) 
by the depolarizing field. Therefore, the lowest subband transition occurs in the dipole left 
FES well case in Fig. 4(b). It causes the abrupt increase of T c /T c0 in Fig. 3. 

Figure 5 shows the effect of the depolarized field strength on the dependence of the ratio 
of ferromagnetic transition temperatures T c /T c0 . We display three different depolarized field 
and apply the bias voltage across the Mn edge-doped and dipole left HDQW. As we expected, 
they have different transition points depending on the F d values. All three lines merge into 
one because the potential profiles of the Mn doped well are not affected by F d , as in Fig. 3. 

In conclusion, we have proposed the DMS/FES hybrid double quantum well structure. By 
using the effects of the spontaneous depolarizing field from the FES well, we can modulate 
the ferromagnetic transition temperature of the DMS well in this system. We calculate 
the Curie ferromagnetic transition temperature in terms of its dependence on the envelope 
functions of carriers at the lowest energy subband. Through the reversal of the depolarizing 
field, we obtain two different ferromagnetic transition temperatures in an applied electric 
field. This result opens the possibility of using ferromagnetic/ferroelectric hybrid quantum 
structures for future multinary spin devices. 
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FIG. 1: (a) Schematic diagram of DMS/FES hybrid double quantum wells. The upper well has 
a 0.5 monolayer of Mn 5-doped layer at the middle of the well (or at the upper edge of the well). 
The lower well is Li-doped ZnO with spontaneous electric polarization, (b) The potential profile of 
the HDQW structures (dotted line) for a dipole up (left) and dipole down (right) cases is shown. 
The confinement potential is Vo = —263 meV with 20% of Mg per unit cell in ZnMgO barriers. 
Vi and V2 are electrostatic potential due to screening charges at the interfaces. 

FIG. 2: (a) Dependence of the ratio of ferromagnetic transition temperatures T c /T c o on bias 
voltage applied across the dipole left HDQW for Mn center-doped (open triangle) and Mn edge- 
doped (closed circle) respectively. The system size is Wi = 10 nm, W2 = 10 nm, B=5 nm. Here 
T c o is T c at F 9 = meV/nm for the Mn center-doped conventional DQW without Li doping, (b) 
Change in the fourth power of the growth direction envelope functions of carriers at the lowest 
energy subband (upper panel), and the potential profile (lower panel) in the HDQW as a function 
of z at F 9 =3 meV/nm (solid) and 7 meV/nm (dashed). 

FIG. 3: Dependence of the ratio of ferromagnetic transition temperatures T c /T c o on bias voltage 
applied across the Mn center-doped HDQW for dipole left (closed circle) and dipole right (open 
triangle) respectively. Inset shows the change of the potential profile due to the reversal of dipole 
polarization in the HDQW as a function of z at F 9 =4 meV/nm. 

FIG. 4: Energy eigenvalues as a function of applied electric fields for hybrid double quantum wells 
for (a) dipole right and (b) dipole left cases. Solid lines indicates the energy upper limit of Mn 
delta-doped well. 

FIG. 5: Dependence of the ratio of ferromagnetic transition temperatures T c /T c o on bias voltage 
applied across the Mn edge-doped and dipole left HDQW for different depolarized electric fields. 



9 



ZnMgO cap 



0.5 ML Mn. 

ZnMgO 




Wi 
B 

W2 




ZnO Buffer 



AI2O3 (0001) substrate 



Kim-Fig. 1 




dipole left at F =4 meV/nm 

dipole right at F g =4 meV/nm 

— i 1 1 1 1 1 1 — 

4 5 6 7 

eV / nm) 



o 
E 

I 

LU 



0) 

E 

LU 



0.00 

-0.05 - 

-0.10-= 

-0.15 

-0.20 - 
'- 

■0.25 - 
-0.30 



-D O O 
=k O o 

-D O O 



O o o 
o o 



(a) 



o o 



° ° « n Q ° ° 
O O « O o o o o o 



° O o 

id o o o 

D ° o o 



o o 

■ 8 8 2° ° ^ 



o o 



o o 



o o 



( 



o 



o o 



o o 







2 



3 



4 



o o 
~ | — r 



Q o 
~ 1 — r 



-9- 



T 

7 



8 



F g (meV / nm) 



0.00 
-0.05 - 

0.10 -i 

-0.15- 
-0.20 



ID 
-0 



-0.25 - D o 



-0.30 



(b) 



§808 



o o 
o o 



o o 
o o 



o o 



o o 



P o 

~D O 



8008 



o 



o o 



800 







T" 

2 



o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


I 

3 




I 

4 




I 

5 



o o 
o o 



o o 
o o 



8 c 



o o 8 o 8 ° [ 



° ? o 



T" 

6 



T 

7 



8 



F g (meV / nm) 



Kim-Fig. 4 



8 



7- 



F =0.001 meV/nm 

d 

F.=0.01 meV/nm 



6- — A — 



F =0.02 meV/nm 

d 



o 

o 



5 - 



4- 



- X 



3- 



2 - 



1 - 



t f 

1 



Kim-Fig. 5 



t ' ♦ t 



F g (meV / ni 




